Landau's quasi-particle mapping: Fermi liquid approach and Luttinger liquid behavior 
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A continuous unitary transformation is introduced which realizes Landau's mapping of the ele- 
mentary excitations (quasi-particles) of an interacting Fermi liquid system to those of the system 
without interaction. The conservation of the number of quasi-particles is important. The transfor- 
mation is performed numerically for a one-dimensional system, i.e. the worst case for a Fermi liquid 
approach. Yet evidence for Luttinger liquid behavior is found. Such an approach may open a route 
to a unified description of Fermi and Luttinger liquids on all energy scales. 
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The behavior of strongly correlated electron systems 
is a fundamental issue with many intriguing but un- 
resolved questions [|]-^. Two antithetic concepts for 
non-symmetry-broken phases govern the discussion: the 
Fermi liquid (FL) [Q and the Luttinger hquid (LL) 
The FL, well-established in three-dimensional (3D) sys- 
tems, is based on Landau's conjecture that the low- 
lying excitations of interacting fermion systems can be 
connected continuously to those of the non-interacting 
Fermi gas: there is a smooth mapping between the quasi- 
particles of the interacting and of the non-interacting sys- 
tem. The LL, well-established in one dimension (ID), 
is based on the observation that collective particle-hole 
excitations exhaust all degrees of freedom. These col- 
lective excitations are bosonic and thus completely dif- 
ferent from the non-interacting fermionic excitations. In 
two dimensions (2D) - relevant for the understanding 
for high-temperature superconductivity - the situ- 
ation is marginal due to the occurrence of logarithmic 
factors. Here we report a non-perturbative prescription 
based on quasi-particle conservation to construct Lan- 
dau's mapping. Evidence is found that LL behavior is 
retrieved with our method. 

FL theory is most successful as phenomenological the- 
ory based on the energy functional 
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where k, q are momenta close to the Fermi surface, 5nk 
the density of quasi-particles - particles above the Fermi 
sea, holes below, e(fc) the kinetic energy and f{k. q) is 
Landau's interaction function. We highlight that Eq. |l| 
impUes the conservation of the number of quasi-particles. 
Microscopic derivations of Eq. |l| are based on perturba- 
tion theory Ordinary electrons or holes are good 
quasi-particles on the Fermi surface k = k-p due to the re- 
stricted scattering phase space. The quasi-particle weight 
< Z < 1 quantifies to which extent, i.e. to which 
fraction, real electrons or holes overlap with true quasi- 
particles. Moreover, the renormalization group (RG) has 
been used in the last decade to give FL theory a firm 



basis Additionally, RG approaches (l^Jl^ are re- 

cently used to understand the 2D Hubbard models in the 
context of high-T^ superconductivity JTsj-p^. All these 
approaches are conceptually based on diagrammatic per- 
turbation theory in the interaction strength. They do 
not employ or construct Landau's mapping except on 
the Fermi surface. Quasi-particle conservation does not 
appear. 

We use continuous unitary transformations (CUT) 
]l6|~p^ to treat interacting Fermi systems 



H = E. 
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J2 r(i^,P, q) ■■ CK+pCK-pCK-q^K+q : (2) 

K,p,q 

where : : stands for normal ordering (in particular 
: Uk '■— Suk) and T{K,p., q) is the vertex function encod- 
ing all scattering processes induced by the interaction. 
We define a CUT which maps the Hamiltonian (||) to 
an effective one HcB which conserves the total number 
of quasi-particles Q := X]fcSigii(|fc| — k-p) : Uk :, i.e. 
[Q, Hcs] = 0. To this end, the flow parameter £ is used to 
denote intermediate Hamiltonians H{£). Starting point 
is H{£ = 0) = iJ; end point is H{£ oo) = H^s. The 
transformation is defined by 



d,H =[rji£),H{£)] 
Vi,j{^) = sign(Qi,i - Qj,j)Hi^j{£) 



(3a) 
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where (3b) with sign(O) — defines the infinitesimal gen- 
erator rj. The matrix elements rii_j and Hij are given 
in an eigen-basis of Q. The flow stops for [Q,H] ~ 
since then 77 = 0. This represents indeed the fix point 
at ^ = 00 as can be shown generally under certain as- 
sumptions The choice ( pb| ) eliminates all parts 
of H changing the number of quasi-particles. The block 
band diagonality in H{£) [19 2Q] is retained which means 



here that the number of quasi-particles is changed only 
by 0, ±2, ±4. Even the 3- or more-particle terms do not 
change it by ±6 or higher. It is very important that (|^) 
is renormalizing since it suppresses matrix elements be- 
tween states with very different energies more strongly 
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than those between similar energies |T£ . Since ap- 
proximations have to be made to solve (| ) this renormal- 
izing property is essential in order not to loose important 
low-energy physics. 

Let us assume that we are able to solve Eqs. ^ exactly 
and that the resulting flow converges. Then the above 
defined CUT represents Landau's mapping of quasi- 
particles. The conservation [Q, Hcs] — allows to de- 
scribe all excitations by a finite number of quasi-particles. 
This is one of the crucial points of the present work. The 
conservation of the number of quasi-particles ensures that 
terms in iJeff consisting of 6, 8, or more fermionic opera- 
tors become relevant only if three or more quasi-particles 
are present in the system. The terms made from 2 and 
4 fermionic operators constitute the 1-particle and the 
2-particle parts of the effective Hamiltonian. They rep- 
resent the dispersion and the effective interaction (cf. Ref. 
Ip^ ) . The basic idea of FL theory is that at low energies 
only the 1-particle and the 2-particle parts are relevant 
even though 3-particle terms are present. Comparing 
Eq. |l| and Eq. | at = oo the vertex function at zero 
momentum transfer {p — q) can be identified with Lan- 
dau's interaction function f{K — p, K + p) — T{K,p,p). 
For non-zero momentum transfer T(K,p, q) can be inter- 
preted as general quasi-particle scattering amplitude. A 
precise one-to-one correspondence, however, to the quan- 
tities known from conventional FL theory Q cannot be 
established since the CUT keeps all interactions local in 
time so that the vertex function F in Eq. |^ has no fre- 





quency dependence. 

So the CUT defined in Eq. || allows to compute Lan- 
dau's interaction function / directly for a microscopic 
model. But there is a caveat. Terms comprising six and 
more fermionic operators do not vanish and they influ- 
ence the 1- and 2-particle terms for < £ < cxd because 
the quasi-particle conservation is achieved only at the 
end of the transformation (£ = oo). The restriction to 
the scattering processes labelled by the vertex function 
F constitutes an approximation. It is related to a 1-loop 
RG calculation for the interaction vertex or, equivalently, 
to the summation of the parquet diagrams [ pT| . Re- 
cently, such RG computations are performed by numer- 
ous groups |l3|- [l5|j2^ for 2D Hubbard models, for which 
results by a CUT ||2^ exist, too. There are three particu- 
larly attractive features of the CUT approach (||). First, 
it is local in the flow parameter £. Second, there is no 
need to take a retardation into account, i.e. no frequency 
dependence occurs because the unitarity of the transfor- 
mation ensures that all intermediate Haniiltonians stay 
hermitian. This is a conceptual advantage over the RG 
procedures used so far since the book-keeping related to 
another dimension, namely time, is avoided. Third, and 
most importantly, a CUT keeps states a all energies. No 
information is lost. This major conceptual advancement 
makes it possible to address also non-universal behavior 
like the full frequency dependence of spectral densities 




FIG. 1. Cuts of r{K,p,q) in Eq. |2| for if = n/4 and p,q e [-n,n]: (a) beginning {£ = 0), (b) middle {£ = 0.5), (c) end of 
transformation {£ = (x>). Black lines stand for processes changing the number of quasi-particles by 0; dark grey for changes by 
±4; light grey for changes by ±2. 



In order to back up the above general results we per- 
formed a numerical calculation for a tight-binding model 
of interacting fermions without spin on a chain with av- 
erage filling n — 1/2. The bare dispersion eo{k) reads 
— cos(fc). The interaction is taken to be a nearest- 
neighbor repulsion implying To{K,p, q) :— V sm[p) sin(g) 
psf . This model is equivalent to an anisotropic, antifer- 
romagnctic Hcisenberg spin chain with 5=1/2 via the 
Jordan- Wigner transformation It is exactly solv- 

able so that much is known about its properties p7| , p8t . 
The system is metallic for U < 1. For larger values the 
discrete translational invariance is broken and a gap is 



opened making the system insulating; here we do not 
consider this phase. The metallic phase itself is a generic 
realization of a LL . The Luttinger parameter can 

be deduced from the comparison of exactly known quanti- 
ties like velocities 31 1 to the predictions of bosonization 
p9| , ^ . So the power-law behavior in the momentum 
distribution at the Fermi edge is known quantitatively 
This power-law behavior and the absence of a 
jump in the momentum distribution is taken as signature 
for the absence of quasi-particles and the prevailance of 
bosonic modes. Thus the ansatz (|^) in terms of fermions 
may seem inappropriate. Yet this is not the case. 
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We use the parametrization in Eq. g with ^-dependent 
energy constant Eq, dispersion e, and vertex function F. 
In the Hmit of an infinite system the starting values are 
T{£ = 0) ^ To, e{l = 0) = -(I + 2l^/7r) cos(fc), and 
Eo{l = 0) = (-(1 + y/2)/7r)A^, where N is the system 
size. The deviations from the bare values result from the 
normal-ordering assumed in Eq. ^. The Eqs. ^ and ^ 
form a closed set of equations once the normal-ordered 
3-particle terms engendered on the r.h.s. of Eq. 3a are 
omitted |2^. We treat the resulting set of functional 
differential equations by discretization, i.e. we solve the 
equations for finite systems. Up to about TV = 50 sites 
could be treated well numerically, which amounts up to a 
set of several thousand differential equations. In Fig. 1, 
cuts of the vertex function F for a given total momen- 
tum K illustrate the transformation at work. Different 
grey scales encode the action of the scattering process 
on the number of quasi-particles. Clearly, the processes 
that change the quasi-particle number are gradually sup- 
pressed (light and dark grey) on £ ^ cxd so that at ^ = oo 
only the processes (black) remain which leave the number 
of quasi-particles unchanged. 
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FIG. 2. (a) Ground state energy per site beyond absolute 
and linear term: exact for N — oo (solid); CUT for TV = 48 
(dashed), (b) Momentum distribution exponent a: exact for 
iV = oo (solid); CUT for A'' = 48 (dashed), details in main 
text; bosonization for TV = oo (dotted). 



In Fig. 2a the ground state energy beyond the abso- 
lute and linear term in V, i.e. the correlation part, is 
depicted. The approximate CUT result agrees well with 
the exact result. Only for V ^ 0.75 a deviation becomes 
discernible; beyond Vsi « 0.94 the approximate CUT 
treatment breaks down because of the omission of the 
3-particle terms. This omission leads to a spurious in- 
stability (subscript si)- Scaling the interaction I4i where 
the spurious breakdown occurs to infinite system size we 
find Vsi{N = oo) « 0.75 which shows that there is a fi- 
nite range of interactions where the CUT restricted to 



the ansatz (g) works. 

More information on the underlying physics is provided 
by the the momentum distribution n(fc) := (uk) which 
can be computed by n{k) = NdE/deo{k). A generic re- 
sult is shown in Fig. 3 (crosses). At first glance, the data 
points towards a fairly large jump. But the points close 
to the Fermi surface bend strongly towards the value 1/2. 
Similar CUT results were obtained for a Tomonaga model 
p2[ . For a quantitative analysis we use a FL fit Q and 
a LL fit as shown in Eqs. ^ and |4^ 

|n(fc)-l/2| « Z/2 -I- A/c (S In(Afc) -t- C) 
|n(fc)-l/2| ^ BAk"" + CAk 



(4a) 
(4b) 

,B,C anda,B,C 



with Afc = |fc— fci?|; the parameters Zk 
are fitted. Checking the applicability of (4b) on finite-size 
data obtained by bosonization for the Luttinger model 
shows that least-square fits are most robust and pro- 
vide good estimates within 5% accuracy for the exponent 
a. In Fig. 3, the excellent agreement of the LL type fit 
and the much poorer agreement of the Fermi type fit are 
salient. This result is not affected by the precise fit pro- 
cedure. 
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FIG. 3. Momentum distributions (data is point-symmetric 
about {tt/2, 1/2)); Fermi liquid fit after Eq. Wat Luttinger liq- 
uid fit with a = 0.076 after Eq. 0. 



By numeric evaluation, we cannot prove the power law 
behavior. For instance, previous results [p^ , which rely 
on additional approximations in the transformation of 
the observable, yielded only a logarithmic correction. We 
find that our numerical data is better described by ( ^ ) 
than by a In-type fit B + a\n{Ak) + CAk just in the same 
way the data in Ref. is better fitted by (^) than by 
the In- type fit. 

In addition, the exponents a in Fig. 2b obtained from 
the fits agree very well with the exact result. Thus the 
excellent agreement of the LL fit in Fig. 3 is not acciden- 
tal. For V < 0.6 the approximated CUT result agrees 
even better with the exact result than does the straight 
application of bosonization. While the bosonization re- 
sult is much too smooth the approximated CUT overes- 
timates the exponent due to the neglect of the 3- and 
more-particle terms. 
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The very good agreements in Figs. 3 and 2b provide 
evidence that the CUT (H,^) targeting on quasi-particle 
conservation captures the relevant LL behavior. This 
comes as a surprise since the concept of holes below and 
particles above the Fermi level as quasi-particles belongs 
in the first place to FL theory. This shows that a power- 
ful RG procedure such as the CUTs allows one to stick to 
the FL type of quasi-particles. An analogous observation 
was made previously for spinous and triplets |Q. It is 
decisive to describe the multi-particle dynamics properly. 

The validity of our mapping onto conserved fermionic 
quasi-particles is strongly supported by recent results 
p5| , ^ for the closely related quantum sine-Gordon 
model. By a similar CUT Kehrein mapped this strongly 
correlated model, which corresponds to an interacting 
fermionic model [^7| , for a wide range of parameters onto 
non-interacting fermions ]38| . 

A description of LL behavior in terms of fermionic 
quasi-particles is intriguing since one can assume that the 
dimensional crossover to higher dimensional FL behavior 
is also captured: for Z? > 1, the quasi-particle description 
should work even better. Existing RG approaches pro- 
vide good insight in the dimensional crossover in terms 
of effective low-energy models [^9|-^. Renormalizing 
CUTs yield a quantitative description on all energy scales 
since no degrees of freedom are integrated out. Besides 
the advantages of the locality in the flow parameter t and 
of the avoidance of retardation, this feature constitutes 
the main progress. In this way, spectral properties with 
full frequency dependence are accessible Q. 

Summarizing, we presented a route based on continu- 
ous unitary transformations (CUT) to the explicit con- 
struction of Landau's mapping of the excitations of the 
interacting system to the ones of the non-interacting 
system. The importance of quasi-particle conservation 
was highlighted. We constructed Landau's mapping nu- 
merically in the worst case, namely for one-dimensional 
fermions. The results indicate that the CUT approach 
captures a finger-print of Luttingcr liquids: power-law 
behavior in the momentum distribution. This finding 
suggests the CUT approach as method to describe the 
crossover from Fermi liquids to Luttingcr liquids on all 
energy scales. 
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